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Introduction 

The  author,  not  being  a  Hydrodynamicist,  was  much  impressed  by  the 
difficulty  experienced  in  finding  procedures  for  computing  or  approximating 
hydrodynamic  forces,  on  bodies,  generated  by  waves.  In  particular,  knowledge 
was  desired  concerning  the  hydrodynamic  forces  and  .moments  acting  on  a  body 
of  revolution  moving  under  waves.  Very  few  papers  were  found  which  attempted 
these  ends.  Some  of  these  are  given  in  the  selected  bibliography.  Most  of 
the  methods  were  valid  for,  or  developed  for,  only  very  restricted  cases  such 
as  balls  or  motion  parallel  to  a  mean  surface,  etc. 

Attention  finally  centered  on  a  paper  by  Cummins  [l3.  However,  the 
exposition  was  found  to  be  lacking,  for  our  purposes,  in  the  following  respects; 

1)  restriction  to  Potential  Flow 

2)  restriction  to  slender  bodies  of  revolution 

3)  restriction  to  motion  parallel  to  the  mean  free 
surface 

h)  restriction  to  constant  velocity  along  an  axis 
of  the  body. 

Vie  shall  consider  these  points  individually. 

The  restriction  to  potential  flow  is  most  easily  treated.  If  anyone 
has  done  something  similar  for  non-potential  flows,  we  have  not  learned,  as  yet, 
of  his  work. 

No  mention  need  be  made  here  concerning  the  limitation  to  treatment 
of  slender  bodies  of  revolution,  except  to  remark  that  with  the  mapping 
techniques,  e.g.  Miles  [2],  available,  this  restriction  is  not  as  great  as 
it  first  appears. 

The  restriction  to  motion  parallel  to  the  mean  surface  is,  we  think, 
removed  by  this  paper.  Ever,  so,  the  orientation  of  the  body  must  remain 
constant.  This  is  not  too  hard  to  live  with  since  this  is  just  what  would  be 
required  for  numerical  integration.  That  is,  if  we  are  integrating,  numerically 
with  respect  to  time,  say,  then  over  a  snail  time  interval  we  would  hold  such 
things  as  orientation  and  velocity  constant. 

The  restriction  to  constant  velocity  along  an  axis  of  the  body  is 
not  needed.  We  shall  still  take  the  velocity  as  constant,  but  will  not  restrict 
it  to  being  along  an  axis. 

The  present  report,  then,  can  not  claim  to  be  anything  more  than  a 
slight  extension  of  Cummin's  results  to  a  situation  which  may  be  useful  in 
computer  applications . 
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The  organization  and  notation  of  this  report  i3  essentially  that 
of  Cummins.  This  should  make  r.v  mparison  easy. 

It  is  recommended  that  this  report  only  be  read  if  the  reader  has 
available  a  copy  of  [l]. 

Coordinate  Systems 

It  will  be  necessary  to  consider  two  systems  of  coordinates. 

The  first  system,  considered  as  stationary  will  be  called  the  earth 
system.  It  will  be  right-handed  and  orthogonal  with  its  origin  at  the  mean 
free  surface  of  the  fluid.  The  positive  direction  of  the  z-axis  is  to  be 
upward  (out  of  the  fluid)  and  is  assumed  parallel  to  the  local  gravity  vector. 

The  fixed  axis  will  be  denoted  by  x,  y,  z. 

The  second  set  of  coordinates  is  also  to  be  right-handed  and  orthogonal. 
This  will  be  called  the  body  system  and  will  be  aliened  in  the  body,  fixed  with 
respect  to  the  body.  The  axis  will  be  denoted  by  x,  y/  z\ 


The  Fluid  and  Its  Boundary  Conditions 

We  consider  a  free  surface  disturbed  by  a  regular  train  of  waves. 

As  a  further  assumption  we  consider  a  linear,  irrotational  wave  theory  for  deep 
water.  Such  a  theory  might,  for  example,  be  specified  by 

V2<p  =  o 

P 

-  =  g(z  +tj)  +  (p ,  +  constant 
P  v 


cp  =  0  on  the  bottom 

t)  =  -  cp,  at  z  =  0  in  which  T)  is  the  surface  elevation. 

o  ^ 

The  linear  irrotational  wave  theory  for  deep  water  is  satisfied  by 
the  velocity  potential  (earth  system) 
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X  =  wavelength 

h  =  waveheight  (crest  to  trough) 
c  =  celerity 
a  «=  phase  angle 

P  ~  direction  of  propagation  as  measured  from  the  x-axis. 
The  fluid  velocities  are  then  given  as 
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For  a  body  moving  toward  the  surface  at  a  constant  velocity  V  we 
obtain  the  following  relation  between  coordinates 
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is  the  matrix  of  direction  cosines  relating  the  systems.  Furthermore 
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In  the  following  we  shall  treat  0._  as  a  constant;  and  for  convenience 

we  set. 


IX  HOKMt  IMVWrTY 
AmifD  PHYSICS  LASOtATOftY 

KV«  MM}  MU1UW 


! '  :■ 


Y  ^  ~ 
c  2n 

6  *  T 


The  potential  function,  with  respect  to  the  moving  system,  is  then 

<P  w  =  <PW  +  •  V  which  may  be  written 

6[%1x'*a32y'ta332'^mt]  ,  , 

%  =  Y  a  cos  6  \_[azix  +aizy  +ai3z 

+UmtjcOS  p  +[  a21x/+a22y  ,+a23Z,+'rrn't]  S^n  P  “Ct  + 

+  (°beP)  *  x 


Hence  we  find  that 
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Singularities 


We  shall  represent  the  body  in  the  same  manner  as  in  [l] . 


is,  we  set 

(li)  H  =  H  i  +  H  j  +  H  k  with 
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Hydrodynamic  Force 


The  Force  is  to  be  computed  via: 


(5)  dF  =  dF,  +  dF^  in  which 

dP;  *  ^  &  >  *,  *' 

dF.  -  -anp  -it-iit-Jv  dy' 

t  v  at  at  3  aT  k>  0X1 

Eqa  [19]  J'rLj?®*0"  used  «»t  of  Cumins  [1],  see  especially  page  5 

We  have  then  the  relation 
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We  find,  from  these,  that 
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(7)  If  3 
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If,  furthermore,  we  consider  the  relations 
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we  may  -write 


dF^  =  -hnp[n.x  (qw)x/+  Hy  (gw)  /  +  Hz  so  that 
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considering  terms  of  first  order  only  (neglect  terms  containing  products 
of  the  variables).  This  may  be  reduced  to 
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Proceeding  similarly  we  find 
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+a22ain  pj  J  a13cos  p+a23sin  p  j  ^  \  \i  iix' 


w  >  dx 


dFZ  14)63  1  y2f13Va33Um  “  COS  p  J  “  [  (anVa21Va3lWm)(a13a3X+a33axi) 

+  2(aX2Va22Va32Wm^ai3a32+a33ai2)  +  Zta^+a^+aaaWj  (2a13a33)  )  U 

/  r  1  r  ‘  J/ 


+  lv2[a23Va33Vm  j  “  [(ailVa2lVm+a31Wm)(a23a31+a33a2l)  +  2(&  12U^a ^ 
■ra32Wm)(a2:,a32+a33a22)  +  2(ai3  Va23Vm+a33Wm)  (  2a2  3a33 )  ^  V  +  ^2  a 33^ 

-  (a13cos  p.a23sin  PXucos  p+Vm  sin  p  .  c)j  -  [(axx^aiVa3i*n)  (a33a3x 
*•  [ax3C0S  P+a23ain  p]  [a71cos  P+a21sin  p]  j  +  (5  ^ 

"  (ai3cos  P'+a23sin  p]  [a12cos  P+a^sin  p]  +  (ai3\+*23\^33\)  U3l 

“  [aX3C0S  P+a23Sin  P]  )  w|  dx' 


33a32 


JM 

l5li 


>  rl 

t'i 

, .  i f< 


,  w 

i  y 


J  vs* 

ill 
8  0 


;is 


j  Pi 

V 
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Let  1^  be  the  x-coordinate  of  the  rear  of  the  body,  and  let  L2  be 
that  of  the  front.  Then 


x 


L2  L2  L. 

\  dF  dx'j  F  »  (  dF  dx':  F  -  I 

J  x  y  J  y  x  J 

L,  1,  L, 


dF  dx". 
z 


Now,  in  each  of  dF  ,  dF  ,  dF  ,  the  only  terras  which  vary  with  x'  are 
x  y  z 

a,  u,  v,  w.  If  we  set  na2  ®  A  then  each  of  dF  ,  dF  ,  dF  is  of  the  fom, 
5  ’  ’  x5  y’  z  ’ 


p6A  (cxu  +  c2v  +  c3w)  dx',  in  which  cx,  c2,  c3  are  constant. 


Hence,  each  of  F  ,  F  ,  F  is  of  the  form 
’  x’  y’  z 

h  fc  , 

P6  (ci  /  Audx'  +  c2  J  Avdx'  +  c3  j  Awdx'  )  . 


This  may  be  written  as 


hncp$ 

X 

We  may  then  write 


,  r2  ez  h  6z 

(c-jcos  p  +  CgSin  P)  |  Ad  sin^dx'-  c  3  J  A  e  cos^dx' 


r.  -  Spnhc 

Px  X 


{  .  (anwmCOS  P+a3iumC0S  P-«3i  ~  2  cos  p  [(au«;a21v;a3lHJ(aua31) 
+  K2um+a22Va32wm^aiia32+a3iaa2)  +  j 

%lVa3  1Vm  "  2[(ailVa21Va3lWm)(a3ia21)  +  (a12\+a22Va32Wm)(a21a32 

L2 

^31^2  )  +  (ai3Um+a23Vm+a33Wm)  (a2ia33+a3la23)  ]  j  sin  P  J  A  e6Z  sin  t  dx  ' 


a13_cos  p+a21sin  p 


J 


w  cos  P+V  sin  p-c  -  (a13u  ->-a21v  +a,,w  )  ( 

m  m  J  AA  m  *si  m  3i  m  \ 


l31 


"[aixcos  P+a2Xsin  p]*)  +  2(a12um+a22v+a32wm)  (a31a32-  [ a^cos  p 
+a21sin  p]  [a12cos  P+a22sin  p])  +  2(a13um+aa3v+a3^m)  (a31a33-[allCos  p 


+a21sin  pj  [a13cos  B+a-^sin  p]  j 


.  6z  , 

A  e  cos  \jr  dx 


/ 
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(22)  Fy  »  ^  |  ^2  a12wmcos  P+a3aucos  p-ca32  - 


COS  pj (aiium+a2XVm+a3lWm) (ai2rt3i 


+a32ail)  +  ^12Um+a22Va32Wm^  ^ai2a32  *  2^13 (ai2a33+a32+ai3>J 
+  (skaV^2^  I"  |(ailVa21Vm+a31Wm)(a22a31+a32a2l)  +  (^V*2^ 


ra  *x  m  m 


+a32Wm^^a22a32^  +  2(ai3Um+a2 3Vm+a33Wm )  (a22a33+a32a23  ) 


J  / 


m  **  ra 

I  ft  6s  .  , 

/  A  e  sin  y  dx 


~  fz  a32V(ai2C°S  P+a22Sin  P)(UmC0S  P+VmSin  ^"C)]  "  [(ailUm  +a21Vm 

+a3iwm)(a32a3i“[ai2C0S  P+a22sin  p]  [axxcos  p+a21sin  p]  )  +  2(aiau+aaavi 
+a32Wm)  (a32  -[aX2C0S  P+a22Sin  P^)  +  (SaSa^ll003  P 

+a22sin  pj  [a13cos  p+a23sin  p  l>])P  cos  \jr  dx7! 

(23)  P_  =  a13wB00s  p*a„^cos  ?-ca33|  -  cos  f 

*^33*11  >  +  2<il12Va»=Va33Wm)(ai3S2'fS3ai3)+2(ai3Va33Va3s''m)(2a13a33)  ] 

*  (2k3V,,»»T»l'f(aiittB4Si%+Siv»H«,3Si«33Si)  * 

'  l\  1  rL= 

+a32Wm^a23a3?  +a33a22  )  +  2  Vm+a33  Wm^  2a23  S33  ^  )  san  P  jAe  Sill  f  dx 

'(2  ^vKa608  P+a23sin  P)(V0S  P'V1"  P-c  “  (axiVa2iVa3iWm)(a33a3i 
-[a13cos  p+a23sin  p]  [a13Lcos  P+a21sin  p])  +  (axau+aaav+a3awm)  (a33a32 
-^3 cos  p+a^sin  p]  [a12cos  p+a22sin  pj)  +  (a^+a^v^^)  (a33 

a  ]\  rL2  gz 

“[*!■». <*0S  P+a23Sin  m  /*.“  sin  \fr  dx'i  . 


A  e  sin  f  cbo  . 


Of  course,  we  must  remember  when  evaluating  the  integrals  that  z,  f,  A 
are  functions  of  x'. 


-  lb  - 
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Hydrodynamic  Moment 

The  analysis  as  done  by  Cummins,  from  Eq,  [IjO]  on  page  7  to  Eq.  [59] 
on  page  10  of  Ref.  [1]  is  carried  over  with  but  little  change.  For  convenience 
we  shall  write  some  of  the  pertinent  lormulae. 


( 2lt )  M  -  +  Mt 

(25)  m£  » x!  ^ x  w + k*py  ^ x  a)± 

i  i 

(26)  Mt  =  p-~  J  $(rxn)do 

s 

To  evaluate  the  surface  integral  we  proceed  as  in  Cummins.  Correspond¬ 
ing  to  his  Eq.  [b8]  we  have 


(27) 


Writing 


dK^  =  pax') 


_d 

dt 


2n 

J  +  tPs) (.j  sin9  -  k  cosG)  d6 


dx 


and 


we  obtain 


<P^  (*',  a  cose,  a  sine)  =  <p'r(x',  0,  0)  -  a(v'cos  6+  w'sin6) 
<P  =  -a(v'cos  6  +  w'sinO)  ,  as  in  Cummins, 


(28)  dKt  =  -2npa2x' 


m  (l/)  j '  ik <v/)  k 


dx'  . 


We  may  write 


imp  (qw  x  p  )  dx.' 


/ 


(29) 


+ 
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This  is 


dM^  =  (  -  x'dF^  z  +  npa  u/w/dx/)j  +  (x'd F^  y  -npa  u'v'dx ')k,  since 

u=  jj  a2u/(w/j-v/k)  and  the  terms  z^dF^  x,  -y,dFg  x  do  not 

appear  because  the  legally  moment  is,  in  this  case,  evaluted  along 
the  x-axis. 


=  x’dF^  z  +  npa2u/w>dx  -  2npa2x7 (w*)  dx^j 
+  (x'dF^y  -  npa2uV/dx/+  2npa2x/  —•  (v')  dx'jk 


We  shall  write  this  as 


dM  =  dM  j  +  dM  k  . 
T  a 


Then  we  find  that 


2  /  /  /  iitt  ~ 

=  npa  <u  w  -  x  (u  w  ,-*•  2v  w  ,+  2w  w  ,+  2an_u,  +  2a0,v, 

L  x  y  z  13  t  23  t 

+  2a33wt  )jdx' 

=  npa2|x'(uv^+  2vV,+  2w'v',+  2a12u,t  +  ?a22vt  +  2a23wt)  -u'v^dx' 


If  as  a  first  approximation  we  take 


uV  = 


u'v'  = 


t  / 

U  W  /  = 
X 


v'w',  = 

y 


-  (a,1Hn+a21Vr,+a31W™)  w#  "  (*,,11  +a  V  ^3  W  ) 

u  m  n  ij  3i  u  m  m  m 

-  (a..u  +a__v  +a,.w  )  v'  -  (a1Pu  +a?Pv  +a~.,w  )  u' 

Ilia  21  ra  31  12  rn  22  m  32  m  ' 

"  ^anum+a2ivm+a3iwm)  V 

'  wy'  >  etc-> 


r*t  Ky»<s  HOnma  u<in>yri 
>*t iHO  P.nSlC5  iAJOIATOat 

in***  kiinwe 


we  find  that 


dM  =  npa2<  x'  (a,_u  +a_.v  +a_.w  )w',+  2(a,,u  +a  v  +a,_w  )  w', 
y  11  m  21  pv  3i  m  x  i2  m  22  m  32  m  y 

« 

+  2(a. ,u  +a,,v  +a-,w  )  w',-  2a, 3u.  -  2a _,v,  -  2a, ,w, 

13  in  23  m  -*3  ra  z  13  t  22  t  32  t 

-  (a„,u  +a,,v  +a_,w  )  u'  -  (a,,u  +a,,v  +a,,w  )  w  >  dx' 
v  13  m  23  rj  33  my  '  n  n  21  in  3i  nr  j 


dM  =  npa2  <  (a,  .u  +a,,v  +a3,w  )  v'  +  (a12u  +a,2v  +a32w  )  u' 
z  1  11  m  21  m  31  m  12  m  22  m  32  m 


i-tou  ‘aoov  »oon  / 

12  m  22  m  32  m 


-  *'  KiVa»iVa3iw,n)  vx'  +  2<aiaV*aaVa«w«)  V 


*  2(ai3Va33Va33Hm)  -  2a33"t  *  2a33Tt  *  2a32»t  f 


A  straightforwaixi,  but  messy,  substitution  for  the  prinied  variables 
will  show  that  each  of  dM^.,  dM^  is  of  the  form 

npa2  (c,  +  c2x')u  +  (c3  +  c4x')v  +  fc5  +  cgx')  w  dx' 

2 

in  which  c2,  c3,  c4,  cs,  cg  are  constants.  On  putting  na  =  A  and 
integrating,  vie  find  that  dKy  and  dMz  have  the  form 

f  h  6z 

p6y  (cx  cos  p  +  c3  sin  p)  A  e  2  sin  ^  dx*  +  (cg  cos  p 


Pz  6z  ht 

+  c4  sin  p)  :  Ax'  e  sin  f  dx  -  /  (c 


/ \  .  - dz  .  ,  /I 
+  c,  x  )  A  e  cos  t  dx 


mt  k>**i  vT'ftiurt 
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.•(4  via*^  VAiu*»o 


The  results  of  the  integration  are 


(35)  My  =  ^  f  -  (ailum+a2ivni+a3i>:m!^ai3c0s  P^3sin  P)  ♦ 

i  At  S2  f 

+a33wrn^aiic0s  P+a2is^nP)j  J  A  e  '  sin  \|r  dx  +  4—  (a^u^+a^v^ 

Lx 

+a3iwm^  ( [ai3a3i+a33aii]  cos  P  +  sin  pj  +  2(a12ujn+a;32vm 

+a32wm)  (  [*13*38 +^33*x^  cos  P  +  [^aSa^^ij  sin  p)  + 

+a33wm)K3a33C0S  P  +  ^a^ssin  P)  -  ^(^cos  p+a^sin  p)  -  2(a33uicos  p 

y  2  2tz 

'ca33+a33Vin  P)  |  Ay/  eT"sin^  dx7  +  (ai1um+a21vin+a3lWm)a33 


+  (*i3  V^a V^aV  «3i  J A  e  ^  cos  *  dx'  +  |L  (auuin+a21vm 
+a3iwm)  (*33%x-  (ai3C0S  P  +a23sin  p]  [ axicos  p  +a21sin  p])+2(al2un 

+a22Vm+a32Wm)  (  a33a32~  [ai3C0S  P+a23Sin  P]  [ ai2C0s  P+a.,2Sin  pj  j 

+  2(a13um+a23Vm+a33Wm)  (a332-  [al3cos  p  +a233in  pf)  + 

+  2a13 cos  p (u^cos  P+vsin  p  -  c)  -  2a33w ^  +  2a23sin  p(u  cos  p 


+v  sin 


i  Ja  2nz  -> 

f  -  c)  i  /  Ax'e  A  cos  ,if  dx'  l 
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vivn  tftwo  iutriwc 


4“  |  (an%+a2iVa3iWm)  (ai2cos  P+a22sin  p)  +  (a2iVa22Vrr 


3j»j  2nz 


+a33Wra^ailC0S  P  +a21sin  P)j  j  A  e 

J  L1 

+a21Vra+a31Wm^  ([ai2a31+a32all]  COS  (?  +  [ 


e  ^  sin  t  <& ' 


4s-  (aniu 


a22a31+a32a 


2]  sin  p) 


+  ^ai2\+a22\+a32\)  (ai2a32cos  P  +a22a32sin  p)  +  2(a13iy*a83vn 

+a33Wm)  (  [ai2a33+a32ai3j  C0S  P  +  [ a22a33+a32a2^  Sin  p )  -  ^(a^COS  p 

L2  2nz 

/  ~X~  / 

Ax'  e  sin  \jr  dx ' 

'  Ln 

,  1  -  L2  2nz 

-  [  (ailVa2xVa3XWm)  a32  +  (a2lVa22Va32\)  ^  /  A  6  X  COS  ^  dx ' 

L1 

+  (a11um+a21vm+a3lwm)  (  a32a31  -  [aX8cos  P+a22sin  p]  [a^cos  p 

+a21sin  p]  j  +  2(a21um+a22vm+a32wm)  ^a32  -  |a12cos  P+a22sin  p  '  j 
+  2(a31um+a3avm+a33wm)  (a32a33  -  [a12cos  p  +a22sin  p]  [a13cos  p+a83sinp]  ) 


^a32  “  !  a12cos  p+a22sin  p  j 


+  2(a  2cos  f  +  a22sin  p)(u  cos  p+v sin  p-c)-2a32w 


2nz 

“X* 

Ax'e  cos  'j.'  <5x 


/ 


■» 


\ 
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